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The electronic structure of noble-gas solids is calculated within density func- 
tional theory's exact-exchange method (EXX) and compared with the results from 
the local-density approximation (LDA) . It is shown that the EXX method does not 
reproduce the fundamental energy gaps as well as has been reported for semicon- 
ductors. However, the EXX-Kohn-Sham energy gaps for these materials reproduce 
about 80 % of the experimental optical gaps. The structural properties of noble-gas 
solids are described by the EXX method as poorly as by the LDA one. This is due 
to missing Van der Waals interactions in both, LDA and EXX functionals. 
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I. INTRODUCTION 



Density Functional Theory (DFT) is nowadays by far the most popular ab initio method to calcu- 
late ground state properties of atoms, molecules and solids. Its generalization to time-dependent 
phenomena'* extended the scope of applied DFT to electronic excitations. DFT's successes rely 
on the construction of accurate and practical approximations to the exchange-correlation part of 
the energy functional. The local-density approximation (LDA) proposed at the very beginning 
of DFT*^^ was and remains a simple tool to use, whose accuracy, however, is surprisingly good 
given the simplicity of the approximation. Over the years, LDA has been applied to practically 
all systems of interest. In most cases, LDA describes ground state properties, such as equilibrium 
structures, adiabatic phonons, to mention only a few, surprisingly well, even for materials with 
strongly inhomogeneous electron densities. On the other hand, attempts to describe energy gaps. 
Eg, or electronic excitations brought rather limited success. 

The case of Eg is a special one because Eg is defined as the difference between the ionization 
potential, /, and electron affinity, A. Hence, Eg can be written in terms of total ground-state 
energies of systems with different numbers of electrons, so it is in principle a ground state property: 

F.g = I -A = E[N+1] + E[N -1]-2E[N]. (1) 

Here E[N+1], E [N] , and E[N—1] are the total ground-state energies of the systems with N + 1, N 
and A'' — 1 electrons respectively. It has been shown^'^ that Eg can be re-expressed in the following 
form; 

Eg = Ef ^ + A,e, (2) 

where 

E^^ = eN+i{N) - CNiN), (3) 

is the Kohn-Sham (KS) gap and A^c is the discontinuity in the exchange-correlation potential 
under adding and subtracting an infinitesimal fraction, oj, of the integer particle number, N: 

Applying LDA to calculate Eg from Eq. (1), one faces the basic problem that LDA docs not 
bind the A'' + 1 atomic electron system. On the other hand, using Eq. (2) within LDA results in 
the equality of the LDA KS gap and the fundamental gap since the discontinuity A^c vanishes in 
this approximation. As is well known, the LDA-KS absolute gap is always much smaller than the 
experimental gap, in some cases, it is even qualitatively wrong. Namely, LDA predicts zero gap, a 
metal, where in nature the system is a semiconductor. 




A most interesting development in density-functional theory during last few years was the appli- 
cation of the exact-exchange KS method (EXX) to the case of crystalline solids. For this method, 
the total-energy functional is given by: 



E[n] = Ts[n] + J nir)Ve,t{r) + \ j j ^yz^T + ^""'"^ + ^""'"^ 



(5) 



where 



Ex[n] = --y I I (6) 
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is the exact exchange energy functional. Ts[n] is the kinetic-energy term of non- interacting elec- 
trons and the correlation-energy term -EcIj^] remains to be approximated, for example, by LDA. 
The orbitals (/)fe in the expression for Ex are KS orbitals, that is the orbitals which minimize Ts [n] 



for a given density, n(r), and are therefore functional of the electron density. For this reason, T5 
and Ex are both impHcit functional of the density, n{r), and the KS equations can be derived: 

2^V2 + Kxt(r) + VH[n]{r) + Vx[n]{r) + Vc[n]{r)j <Pka{r) = Eka<f>ka{r), (7) 

where the KS exchange potential = 6Ex[n]/Sn{r) is obtained by the chain-rule differen- 

tiation: 
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Here Vifs[n](r) is the total KS potential corresponding to the density, n{r). Neglecting correla- 
tion, the total-energy functional has apparently the form of a Hartrec-Fock expression for the total 
energy. However, the EXX method and the Hartree-Fock method differ because the Hartree-Fock 
and KS orbitals (j>k are not the same. In the Hartree-Fock method these orbitals obey an equation 
with a non-local potential (the Fock operator), whereas in the KS method they are determined by 
a KS Hamiltonian with a local, multiplicative potential. 

For more than two decades, the EXX method, sometimes under the name optimized effective 
potential (OEP), has been applied to atoms and molecules. ^'^ Total ground-state energies were 
always found to be extremely close to the Hartree-Fock values while the single-particle spectrum 
- yielding, for example, a Rydberg series for finite systems - is strongly improved over Hartree- 
Fock. Only recently, however, has this method been applied to bulk semiconductors, insulators 
and metals. A surprising result of EXX calculations for sp-semiconductors is that the KS 
energy gaps between occupied and unoccupied states are very close to experimental gaps. This 
contrasts typical LDA results in which the gaps are always too small. In the EXX case, a larger 
gap than in LDA is to be expected because the EXX potential is self-interaction free and, thus, 
binds more strongly than the LDA potential. For example, the EXX-KS electron binding energies 
in atoms are greater in magnitude than their LDA counterparts and closer to experiment. Since 
the occupied valence states generally have a greater self-interaction contribution than the more 
delocalized unoccupied states, the EXX-KS gap should widen as compared to the LDA. If the 
fundamental gap is evaluated from Eq. (1) with total energies in Hartree-Fock approximation, one 
obtains 

Ef^ = er+,{N)-er{N). (9) 

Evaluating, on the other hand, the gap in EXX, Eqs. (2) and (3) yield 

E™ = -^^(iV) - e^^^-^^(Ar) + A. (10) 

where A^; is the X-only discontinuity. Assuming that, like for finite systems, the total energies in 
Hartree-Fock and the total energies in EXX are very close, one would expect that Eqs. (9) and 
(10) yield very similar values for the band gap. This is indeed the case.^^ What is surprising and 
still not well understood is why the EXX-KS gaps alone, i.e. without adding the discontinuity, are 
so close to the experimental gaps for sp-semiconductors. The prevailing belief is that the exact 
KS gaps are smaller than the true gaps for solids. As there are no exact KS potentials for solids 
available, this belief is yet unverified. It is not known how large the discontinuity, A,,.c, given by 
Eq. (4), of the exact exchange-correlation density- functional is. In any case, the EXX method 
is the only KS method so far yielding the electronic structure close to experiment, provided one 
does not add the discontinuity which should, in principle, be considered. This success has some 
important practical consequences. One of them is the possibility of using the one-electron energy 
spectrum as a first, approximate description of excitations and as a fast tool to interpret experi- 
ments. Another, perhaps more important, is a better starting point for more accurate calculations 
of excitations via the time-dependent density-functional theory. 

In this paper, we explore the performance of the EXX method for materials very different from 
sp-semiconductors, namely, the crystalline noble gases. Solid Ne, Ar, Kr and Xe are special sys- 
tems since these materials are composed of almost independent atoms. In fact, the shape of the 



charge density of superimposed isolated atoms is rather similar to the charge density of the solid 
and wc might expect the energetics to be similar as well. The individual tightly bound close-shell 
atoms bind very weakly with one another. Therefore these systems are a diSicult case for LDA. 
Indeed, as will be seen in Table IV, the LDA description of the equilibrium structural properties 
for these solids is much less satisfactory than usual. These solids are large-gap insulators with their 
experimental energy gaps ranging between 21.4 eV for neon and 9.8 eV for xenon. Furthermore, 
they have a large exciton binding energy on the order of several eV. Since the noble-gas solids are 
a loosely bound ensemble of atoms, a comparison between various electronic properties in the solid 
phase and in isolated atoms is meaningful and interesting. We will investigate both cases within 
the LDA and EXX methods. As will be seen later, in the atomic case the EXX without correlation 
resembles highly accurate Kohn-Sham results when they are available. Adding LDA correlation 
to the exact exchange potential usually worsens the quality of results. However, for the physical 
quantities involving total-energies differences an account of correlation seems to be necessary. 

Our results show that the EXX approximation, with and without LDA correlation, widens the 
LDA Kohn-Sham gaps of noble-gas solids by about 1.3 - 3.4 eV. However, the KS energy gaps 
are not as close to the experimental gaps as reported for sp-semiconductors. Along with previous 
observations that the EXX-KS gaps for diamond^'^ and MgO^ are smaller than experimental gaps, 
our results suggest that, at least for large-gap insulators, the EXX theory does not perform as well 
as it does for sp-semiconductors. 

It has been suggested that for atoms the exact KS gap between the highest occupied and the 
lowest unoccupied state represents an approximation of the optical gap rather than quasiparticle 
gap^®'^^. The latter represents a nonneutral excitation or a total-energy difference between the N 
and ± 1 particle systems while the former is the gap between the ground and excited states of 
the A'^-electron system. This conjecture has been explicitly verified in the cases where almost exact 
Kohn-Sham atomic potentials have been calculated^*''^^. Our EXX results for noble-gas atoms 
(Section II) support this conjecture as well. In regard to the solids, the picture is not as clear. For 
semiconductors, the exciton binding energy is very small, of the order of meV, the observed agree- 
ment between the EXX-KS gaps and experimental gaps could support the conjecture of Refs. [ 
16, 17] as well. However, unless better approximations to the correlation-energy functional are 
known this remains as a rather inconclusive speculation. Our results for noble-gas solids, systems 
where optical and quasiparticle gaps differ by several eV and where electronic structure in the 
solid phase is similar to the structure of energy levels in isolated atoms, could shed some light on 
the validity of this conjecture. What we have found is that the EXX gaps for noble-gas solids are 
- unlike the case of semiconductors - significantly smaller than the quasi-particle gaps. However, 
they are also smaller than the experimental optical gaps, reproducing about 80 % of their values. 
A fundamental question remains then: would an exact treatment of correlation align the KS and 
optical gaps for noble-gas solids as well, or does correlation in the solid phase have a qualitatively 
different character and role than for finite, atomic systems, so that the conjecture of Refs. [ 16, 17] 
is not correct for solids ? 

The next section is devoted to results for isolated atoms. Electronic properties of the noble-gas 
solids are presented and discussed in Section HI. 

II. NOBLE-GAS ATOMS: NE, AR, KR AND XE 

Since noble-gas solids consist of loosely boimd individual atoms, it is instructive to start the 
investigation by analyzing the properties of isolated atoms. In addition, the density-functional 
methods for atoms are more advanced than in extended systems and can be used as a test for var- 
ious approximations. This is for two reasons. First, highly accurate KS potentials are available for 
a few light atoms like Ne^*'^^ and Ar^^. Second, for finite systems it is computationally feasible to 
calculate the fundamental gap, Eq. (1), from total-energy differences. Note that for finite systems, 
the fundamental gap as given by Eq. (1) is usually called the chemical hardness. For simplicity, we 
keep the term "fundamental gap" for both finite and infinite systems. In the following we present 
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several results for the noble-gas atoms and draw some conclusions about the validity of certain 
approximations. 

In Tables I and II we present the KS eigenvalues for the Ne and the Ar atom, respectively. 
Results using LDA, pure EXX (EXX), EXX plus LDA correlation (EXXc) and highly accurate 
{almost exact) KS potentials for Ne and Ar atoms are presented. The column denoted QMC in Ta- 
ble I shows eigenvalues resulting from the KS potential obtained by Umrigar and Gonzc"'^* through 
a Quantum Monte Carlo calculation of the ground-state density of the Ne atom. The column 
denoted CI in Tables I and II shows the results obtained by Morrison and Zhao^*^ from highly 
accurate densities which were calculated through the configuration-interaction approach. The col- 
umn called Expt. contains the negative of the experimental first ionization potential of each atom 
which should be compared to the uppermost occupied KS eigenvalue. The column also contains the 
experimental binding energies of a few excited states. The latter are defined as the negative of the 
difference between the experimental ionization potential and the experimental optical energy gap 
between the ground state and the relevant excited state. For the optical gap we take a weighted 
average of singlet and triplet gaps. As for the occupied levels, we report only the experimental 
ionization potential, i.e. the binding energy of the highest occupied state. This is because it is 
unclear how the lower KS eigenvalues compare with the physical binding energies. However, for 
unoccupied states, it has been suggested that the KS eigenvalues are a good approximation to the 
experimental binding energies of the excited atoms^^'^''. This idea has been explicitly verified for 
helium^^, beryllium^^ and neon atoms^^; and some plausible arguments have been put forward to 
explain why.^^ 

Tables I and II show that the LDA KS eigenvalues are smaller in magnitude than their EXX and 
EXXc counterparts. The magnitude is largely a consequence of the self-interaction error inherent 
in the LDA method. Another consequence of this self-interaction error in LDA is the absence of 
none higher but the first bound, unoccupied state. These facts are well known. What is more 
interesting for this study is that EXX eigenvalues for Ne and Ar are deeper than the exact ones. 
These facts are well known. What is more interesting for us is that EXX eigenvalues for Ne and Ar 
are deeper than the exact ones. This means that the EXX method binds atomic states of Ne and 
Ar too strongly and exact correlation would correct for this over-binding. When LDA correlation 
are added to the exact exchange (EXXc method) the situation worsens. For example, the EXX 
uppermost occupied eigenvalues are roughly 1.5, 0.3 eV deeper than the exact eigenvalues for Ne 
and Ar respectively. In the EXXc method the difference grows to 3.2 and 1.7 eV respectively. It is 
interesting to note that for occupied states this expected behavior of exact correlation of the KS 
theory is similar to the role of correlation in the many-body Green's functions approach. There, 
the many-body correlation shift the eigenvalues, or rather the poles of the one-particle Green's 
functions, in the opposite direction as the Hartree-Fock shift. The well known effect is that bind- 
ing energies of occupied states are much too strong within Hartree-Fock. It turns out that the 
EXX KS occupied states are also more bound than the expected exact KS states, but the effect is 
much smaller than in the Hartree-Fock theory. The analogy brakes down for empty states. In this 
case, EXX results in still deeper unoccupied states than in experiment and adding LDA correla- 
tion makes them even deeper. On the other hand, Hartree-Fock theory hardly binds at all. This 
striking difference between the EXX and Hartree-Fock theories is because Hartree-Fock only has 
an asymptotic — 1/r potential for the occupied states and has an exponentially decaying potential 
for the empty states. In contrast, the KS Vxc potential has a Coulombic tail for all states. One 
says, that the Hartree-Fock potential is self-interaction free only for occupied states, whereas the 
KS potential (exact and EXX) is self-interaction free for all states. 

For the heavier atoms Kr and Xe, we perform a full relativistic OEP calculation as relativistic 
effects should be important. In these cases, we do not have exact KS potentials or eigenvalues with 
which to compare the results. However, a comparison with the experimental first ionization poten- 
tial shows that, when the spin-orbit interaction is neglected, the exact exchange calculation gives a 
slightly over-bound uppermost occupied level. Here again, LDA correlation lower the eigenvalues 
too far. It is clear that for all noble-gas atoms adding LDA correlation to exact KS exchange 
deteriorates the one-electron properties. 



So far we have compared LDA, EXX and EXXc KS eigenvalues with the exact ones, where avail- 
able, and with the experimental first ionization potential and binding energies of unoccupied states. 
For finite systems, it is possible to relate certain excitations with total-energy differences. This is 
the case for the first ionization potential / = E[N — 1] ~E[N], electron affinity A — E[N] —E[N+l\, 
and the energy gap Eg, Eq. (1). The electron affinity, A, as defined by a total energy difference is 
zero within LDA theory for most atoms. This is because in LDA the corresponding A'^ + 1 electron 
system is not bound. In the case of noble-gas atoms, however, the experimental affinity does van- 
ish, and the LDA result is fortuitously correct. Since the affinity vanishes, the gap must equal the 
ionization energy. The same argument is valid for EXX. On the other hand, EXXc gives a small 
but finite value for the affinity. Table III shows the atomic energy gaps Eg, both experimental 
and calculated from total-energy differences within LDA, EXX and EXXc. These energy gaps 
are compared to the KS eigenvalue differences eN+i{N) — €n{N). In addition, the (triplet-singlet 
averaged) first optical gap is presented Aopt, as well as a calculated in each method total-energy 
optical gap A* = E[N]* - E[N]. Here. E[N] is the total ground-state energy and E[N]* is the 
self-consistent total energy of an excited state in which there is a hole in the uppermost p-shell 
and an electron in the next s-shell. 

First, we note that the KS gap rather poorly approximates the experimental fundamental gap. 
The average deviation over the four elements between the KS and true gaps is 5.9, 4.0, 3.1 eV 
for LDA, EXX and EXXc respectively. As shown by the data, LDA KS gaps deviate the most. 
When expressed in porccnts. LDA. EXX and EXXc KS gaps account for 63 %, 75 % and 80 % of 
the experimental fundamental atomic gaps respectively. For Ne (Table I), the highly accurate KS 
energy gap amounts to 77 % of the experimental one. For Ar (Table II) we could only estimate 
the exact KS gap to be of about 68 % of the experimental one. The situation changes when we 
compare KS gaps to experimental optical gaps. It turns out that LDA, EXX and EXXc reproduce 
the atomic optical gaps with 13 %, 3 % and 9 % accuracy respectively. The highly accurate KS 
gap of Ne amounts to 99.9 % of the optical gap. In particular, the good agreement of the EXX 
method with experiment is to be noted. An interesting question is whether similar trends take 
place in the solid phase. We might expect the answer to be yes if the solid is composed of a bunch 
of weakly interacting atoms. 

When we compare the experimental fundamental gap of noble-gas atoms with the calculated 
total-energy differences (Eq. 1) the agreement is fairly good. Across four elements, LDA, EXX 
and EXXc reproduce the experimental atomic gaps on average to 4 %, 6 % and 3 % accuracy re- 
spectively. Performing a similar although less rigorously justified total-energy calculation of optical 
gaps the agreement amoimts to 5 %, 5 % and 4 % for LDA, EXX and EXXc respectively. What 
is remarkable is that for a total-energy difference calculations accounting for correlation seems to 
be very important. In fact, the best total-energy difference results come from the EXXc method, 
the worse with the pure EXX one. 

In the next section we will investigate whether similar trends take place in the solid phase. 



III. NOBLE-GAS SOLIDS: NE, AR, KR AND XE 

We performed our calculations for noble-gas solids within the pseudopotential and plane-wave 
formalism. For each approximate E^c (LDA, EXX and EXXc), a pseudopotential was generated 
using the same functional^^ . We followed the EXX plane- wave formalism developed by Gorling^^ 
and Stadclc et al.^"^. This formalism had been applied in our previous publication^'''. The plane- 
wave cutoff ranged between 50 and 120 Rydbergs depending on the material and whether the 
structural or electronic properties were investigated. For solid krypton and xenon, the spin-orbit 
interaction was taken into account perturbatively. 

In Table IV the experimental and calculated equilibrium fee lattice constants are presented. The 
percent deviations from experiment are also given. As already noted, noble-gas solids resemble 



loosely bound isolated atoms. For such systems the total energy only weakly depends on the in- 
teratomic distance. The energy-volume curve is very flat and the system hardly binds. In the 
absence of stronger interatomic interactions, a theoretical determination of equilibrium properties 
is subject to a rather large uncertainty. Usual convergence criteria for total-energy determination 
can easily become insufficient and minor computational details, like for example, details of the 
pseudo-potential construction can matter. For all these reasons we have checked our LDA results 
against independent results from publicly available LDA codes. ^^'^^ 

As shown in Table IV, the structural equilibrium properties of noble-gas solids are very poorly 
described by the all methods. For example, the error in the estimate of the lattice constant for 
neon is more than 13 %. This discrepancy is unusually large; for most solids, LDA gives a lattice- 
constant estimate within about 1-2 %. For Ar, Kr and Xe, the LDA discrepancy is on the order of 
5 %. Although not as dramatic as for Ne, this disagreement is still much larger than usual. Using 
the EXX method without correlation does not help much; the disagreement with experiment is 
just as large as in the other two cases. In contrast to the LDA which underestimates the lattice 
spacing, the EXX method overestimates it for Ar, Kr and Xe. With the exception of neon, the 
combined method, exact-exchange plus LDA correlation (EXXc), is the closest to experiment, but 
it is still not perfect. This result is similar to what was seen in the previous section's atomic 
calculations. For total energy differences, the EXXc method gives the most accurate results. It 
is worth noting, that none of these approximations properly account for the long-range Van der 
Waals interactions that are responsible for the binding of noble-gas solids. It is not surprising that 
structural properties differ so much from experiment. 

Figures I-IV show the band structures along the L-F-X directions for Ne, Ar, Kr and Xe with Kr 
and Xe being treated rclativistically. The solid lines represent the EXX band structure, the dashed 
lines correspond to the LDA one. In Tables V-VIII the KS energies at high symmetry points F, X 
and L are presented. In Table IX the KS energy gaps are compared with the experimental absolute 
energy gaps and with the corresponding optical gaps for solid Ne, Ar, Kr and Xe. It is important 
to keep in mind that for this special class of elements, the description of the electronic structure 
provided by LDA, EXX and EXXc might be expected to work less accurately in the solid phase 
than for isolated atoms. This is because in addition to only approximating short-range correlation, 
our functionals do not account for the long-range Van der Waals effects. 

Nevertheless, it is interesting to consider how well LDA, EXX and EXXc describe the electronic 

structure of noble-gas solids. In LDA. the KS gaps are on average 55 % of the experimental gaps. 
This is a typical result for solids, and though slightly less accurate than for the atomic limit given 
by LDA. For noble-gas solids EXX and EXXc reproduce about 68 % of the fundamental gap, 
whereas in the atomic limit the KS gap ranged on average to 74 % and 79 % of the experimental 
fundamental gap for EXX and EXXc respectively. The fact that the EXX and EXXc KS gaps in 
noble-gas solids are not very close to the experimental fundamental gaps is a central result of our 
investigation. The atomic results suggested that the EXX-KS gaps should be close to the optical 
gaps. Since the noble-gas solids have a large exciton binding energy, the optical gap is appreciably 
smaller than the fundamental gap. What we observe in Table IX is that the EXX and EXXc KS 
gaps are still significantly smaller than optical gaps in noble-gas solids. They amount to about 80 
% - 81 % of the experimental al-gaps, whereas the agreement is 97 % in the atomic limit. 

The valence bands show very little dispersion in Figs. 1 through 4. Furthermore, the energy dis- 
tance between occupied s and p bands is very close to the KS energy gap between s and p valence 
states in the isolated atoms. This confirms the popular picture that rare-gas solids are composed 
of almost undisturbed atoms. However, the conduction bands exhibit a rather pronounced disper- 
sion. It would be interesting to verify this result experimentally. Clearly, an indirect confirmation 
is the fact that the experimental optical gap of the solid is always larger than the optical gap in 
the corresponding atom: The hole and the electron attract each other strongly if they are both 
localized on the same atom. In the solid, on the other hand, the hole is similarly localized as 
in the atom, while the electron in the conduction band is more delocalized, leading to a weaker 
interaction with the hole. 
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Another remarkable result is that the total valence-band width resulting from LDA is smaller 

than the one from EXX for neon. For argon, both widths are comparable. For krypton and xenon, 
the EXX width is smaller. We recall that previous EXX calculations have shown that for sp- 
semiconductors, EXX leads to narrower total valence-band widths than in LDA. Only for diamond 
was the opposite observed. 



IV. CONCLUSIONS 



We have applied the exact exchange method within density-functional theory to the noble-gas 
solids, Ne, Ar, Kr and Xe. It was previously shown for He, Be and Ne atoms that the Kohn-Sham 
energy gap coming from a nearly exact KS potential is an excellent approximation to the atomic 
optical gap but not to the fundamental (quasiparticle) gap. The EXX-KS gaps for these atoms are 
also in very good agreement (3 % in average) with experimental optical gaps. A central question of 
our investigation was whether the same holds for the noble-gas solids. It turns out that, in contrast 
to previous results for sp-semiconductors, the EXX-KS gaps in noble-gas solids are appreciably 
smaller than the experimental fundamental gaps. Moreover, they are also smaller by 20% than 
the experimental optical gaps. The results of our investigation clearly show that the EXX method 
does not provide a KS band structure that agrees equally well with experiment for semiconductors 
and insulators. 
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TABLE I. Neon-atom energy levels (in eV). The Kohn-Sham energies are from LDA, exact exchange 
(EXX) and exact exchange with LDA correlation (EXXc). Column QMC gives the eigenvalues obtained 
with the almost exact Kohn-Sham potential of Ref. [18]. Column CI presents results of almost exact 
Kohn-Sham calculation of Ref. [20] . 
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TABLE II. Argon-atom energy levels (in eV). The Kohn-Sham energies are from LDA, exact exchange 
(EXX) and exact exchange with LDA correlation Column CI presents results of almost exact Kohn-Sham 
calculation of Ref. [20] . 
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TABLE III. Fundamental energy gaps Eg=I-A and optical gaps from experiment and calculations in 
neutral atoms Ne, Ar, Kr and Xe. Aopt is the multiplet-averaged experimental transition energy from the 

ground state to p''s^ state. A* is the calculated total-energy difference between the excited atom in the 
{p's^) configuration and the ground state. lS.tKohn-sham is the Kohn-Sham gap. 
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TABLE IV. Equilibrium cubic lattice spacing (in a.u.) from experiment and calculations. Percents 
show the diviations of the lattice constant from experiment. 
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TABLE V. 


Nc-sohd Kohri-Sham eigenvalues 


in eV from LDA, EXX and EXX plus 


LDA correlation 


(EXXc) at high symmetry points. 
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TABLE VI. Ar-soUd Kohn-Sham eigenvalues in eV from LDA, EXX and EXX plus LDA correlation 
(EXXc) at high symmetry points. 
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TABLE VII. 


Kr-solid Kohn-Sham eigenvalues in eV from LDA, EXX and EXX plus 


LDA correlation 


(EXXc) at high symmetry points. Spin-orbit sphttings included. 
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TABLE VIII. Xe-soUd Kohn-Sham eigenvalues in eV from LDA, EXX and EXX plus LDA correlation 
(EXXc) at high symmetry points. Spin-orbit splittings included. 
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11.32 
8.16 
6.47 
5.26 
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6.69 



14.76 
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11.6 
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12.2 
10.2 
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TABLE IX. Calculated and measured energy gaps in noble-gas solids in eV. E^°^, Ef and Ef 
are Kohn-Sham gaps from LDA, pure EXX and EXX plus LDA correlation respectively. is the 

experimental fundamental gap. A is the experimental optical gap. 



1 Q 



FIGURES 



Fig. 1. Band structure of Ne along L-F-X directions calculated within EXX 

(solid lines) and LDA (dashed lines). 
Fig. 2. Band structure of Ar along L-F-X directions calculated within EXX 

(solid lines) and LDA (dashed lines). 
Fig. 3. Band structure of Kr along L-F-X directions calculated within EXX 

(solid lines) and LDA (dashed lines). Spin-orbit splitting included. 
Fig. 4. Band structure of Xe along L-F-X directions calculated within EXX 

(solid lines) and LDA (dashed lines). Spin-orbit splitting included. 
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